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Combining the equivalence between cosmological particle creation and an effective viscous fluid 
pressure with the fact that the latter represents a dynamical degree of freedom within the second- 
order Israel-Stewart theory for imperfect fluids, we reconsider the possibility of accelerated expansion 
in fluid cosmology. We find an inherent self-limitation for the magnitude of an effective bulk pressure 
which is due to adiabatic (isentropic) particle production. For a production rate which depends 
O^ ' quadratically on the Hubble rate we confirm the existence of solutions which describe a smooth 

^\ I transition from inflationary to noninflationary behavior and discuss their interpretation within the 

0^ ' model of a decaying vacuum energy density. An alternative formulation of the effective imperfect 

fluid dynamics in terms of a minimally coupled scalar field is given. The corresponding potential is 
discussed and an entropy equivalent for the scalar field is found. 
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I. INTRODUCTION 



It is the essential feature of the second-order thermodynamical theories of Miiller M, Israel and Stewart [0-0, 
. _ and Pavon et al. [g|,|| that the bulk and shear viscous pressures and the heat flux, which describe the deviations 
J~. . of a relativistic fluid from local equilibrium, become dynamical degrees of freedom on their own and obey "causal" 
QQ ■ evolution equations guaranteeing subluminal propagation speeds of thermal and viscous perturbations. By now, it 
tH- ■ is generally agreed that these second-order theories should replace the traditional first-order theories by Eckart and 
f^ ' Landau and Lifshitz, which suffer from serious drawbacks concerning causality and stability. 

In homogeneous and isotropic models of the universe a bulk viscous pressure is the only possible dissipativc phe- 
CTN nomenon. The role of bulk pressures in cosmic media is twofold. Firstly, they generally occur if different components 
^J I of the cosmic substratum are coupled. Due to their different internal equations of state the cooling of the subfluids 
r^ . with the expansion of the universe is different, resulting in a tendency of the system to move away from equilibrium. 
Oh These different cooling rates of the components give rise to a bulk viscous pressure of the cosmic medium as a whole 
O • [0nlll ■ Secondly, a nonvanishing bulk pressure of the cosmic fluid may be the consequence of particle number non- 
Vh ■ conserving interactions. In particular, one may think here of quantum particle production out of the gravitational 
c/j [ field O-llJ] . The manifestation of cosmological particle creation as an effective bulk pressure has initiated a series of 
Ci • papers discussing a fluid phcnomenological approach to cosmological particle production ||lq-p3[ . The present paper 
will mainly explore this second aspect as well. 

Special attention has been devoted to "isentropic" (or "adiabatic") particle production ||l7| , |l8|| , i.e., production of 

perfect fluid particles. Here "isentropic" means constant entropy per particle. There is, however, entropy production 

^ , due to the enlargement of the phase space of the system since the number of perfect fluid particles increases. The 

■ - ' condition of isentropic particle production establishes a simple relationship between the particle production rate and 

the viscous pressure which in most application up to now has been investigated in connection with an Eckart-type 

expression for the entropy flow vector (see, however, ||23| , |34| ) . 

Because of its generally negative sign many authors have addressed the question whether a cosmological bulk pres- 
sure may become sufficiently large to induce a phase of accelerated expansion ( ||35|- pd| , p3| , p5| ) . The basic shortcoming 
of most previous attempts to establish scenarios of bulk-viscosity-driven inflation for cosmic media with conserved 
particle number is an apparently unphysical behavior of the thermodynamic quantities. The point is that the non- 
equilibrium viscous pressure has to be of the order of or even to exceed the equilibrium pressure. This assumption 
is in itself problematic and difficult to justify. Moreover, as a result of the back reaction of a viscous pressure of 
this order on the temperature, the universe heats up while at the same time the cosmic matter is diluted [ p5[ . It is 
difficult to imagine interaction processes with these consequences in particular in the extreme case of an exponential 
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expansion [ pOp^pj] ]. In this respect the causal second-order theory is not really superior to the traditional Eckart 
theory. 

The situation is different, however, if the fluid particle number is allowed to change, i.e., if cosmological particle 
production processes are taken into account. The particle production rate may be regarded as an additional parameter 
for the fluid dynamics. Under such conditions the corresponding effective bulk pressure is not necessarily small. For 
example, the production rate may be fixed by symmetry requirements, corresponding to specific self-interactions of 
the cosmic medium | |3l| , |3^ |. 

In the present paper we follow a different strategy. We combine cosmological particle production with an entropy- 
flow expression of the Israel-Stewart type. However, our considerations continue to rely on the concept of isentropic 
particle production and therefore retain the simple relationship between creation rate and viscous pressure, mentioned 
above. This procedure allows us to replace the bulk pressure in the second-order evolution equation by the corre- 
sponding particle production rate and to make the latter quantity an independent degree of freedom. The particle 
production rate is dynamized and no longer a free parameter. This provides the framework in which we reconsider 
the question whether phases of inflationary expansion may be consistent with a fluid cosmological description. It 
turns out that this setting avoids the drawbacks of "bulk viscous inflation" without particle production. We establish 
a model where the universe starts in a de Sitter phase with finite, stationary values of density and temperature, and 
subsequently evolves smoothly to a standard FLRW stage. We point out similarities and differences to the scenario by 
Gunzig et al. 129], which is based on the assumption of a deca ying cosmological "constant". In particular, within the 
causal, second-order theory the corresponding "deflationary" ^^ solution is subject to additional thermodynamical 
constraints which represent a self-limitation concerning the admissible magnitude of the effective viscous pressure. 

Moreover, we clarify how the effective imperfect fluid picture is related to the dynamics of a scalar field. While 
the scalar-field-based inflationary dynamics relies on a "slow roll" period with "adiabatic supercooling" , followed by 
a "reheating" phase during which all the entropy is produced, the fluid picture implies entropy production already 
during the de Sitter phase. The interrelation between the fluid and the scalar field pictures will allow us to express 
the growing entropy (due to particle production) in a comoving volume also in terms of the scalar field variable. 

The paper is organized as follows. In Sec. II we summarize (and somewhat generalize) the second-order theory of a 
homogeneous and isotropic, spatially flat bulk viscous universe and point out the problems inherent in "bulk viscous 
inflation" without particle production. Section III is devoted to the dynamics of a universe with "isentropic" particle 
production on the basis of a second-order expression for the entropy flow vector. The creation rate is regarded as a 
dynamical degree of freedom. We consider linear and quadratic dependences of the particle production rate on the 
Hubble parameter and discuss the self-restriction of this rate as a consequence of the positivity requirement for the 
corresponding relaxation time. For the case of a quadratic dependence an inflationary model of the universe with 
a smooth transition to a Fricdmann-Lemaitre-Robertson- Walker (FLRW) stage is established and its relation to a 
scenario with decaying cosmological term is discussed. The relation between fluid cosmology with isentropic particle 
production and an alternative scalar field description is clarified in Sec. IV. In Sec. V we summarize our results. 
Units have been chosen so that c = ks = h = 1. 

II. THE VISCOUS UNIVERSE 

We assume the cosmic matter to be characterized by the energy momentum tensor of a bulk viscous fluid, 

T''' =pu'u'' + ip + n)h''' , (1) 

and the particle flow vector 

N' = nu' . (2) 

Here, p is the energy density, p is the equilibrium pressure, n is the particle number density, w* is the fluid four- velocity 
in the Eckart frame, and h'^'^ = g'^'^ + u^u^ is the spatial projection tensor. The quantity tt denotes that part of the 
scalar pressure which is connected with entropy production. The conservation laws N""^ ~ and T"^ = then imply 

n -(- en = (3) 

and 

p = -Q{p + p + ^) , (4) 

respectively. Taking into account second-order deviations from equilibrium, the entropy flow vector S*" takes the form 
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where s is the entropy per particle, T is the fluid temperature, and t is a relaxation time. The coefficient of bulk 
viscosity is denoted by C,. For the entropy production density one obtains 



The simplest way to guarantee 3°"^ > is a generalized linear relation 
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which implies the evolution (see, e.g., ]25[| and references therein) equation 
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for TT, as well as the expression 
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for the entropy production density. It is the basic feature of the causal non-equilibrium theory that the viscous 
pressure tt becomes a dynamical degree of freedom on its own and is subject to an evolution equation (ph . For r — > 
Eq.(||) reduces to the corresponding algebraic relation tte = — C© of the Eckart theory (subscript E). 
We will assume the equilibrium variables p and p to obey equations of state 



p^p{n,T) , p^ p{n,T) . 

The temperature behaviour of a bulk viscous fluid is governed by ]l8|,p^,p5p] 
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For TT — the temperature law (nl]) reproduces the well known T oc a~^ behaviour in a radiation dominated FLRW 
universe, while we recover T oc a^^in the matter dominated case. For a viscous fluid the behaviour of the temperature 
depends on tt. Since tt is expected to be negative, the second term in the bracket on the right-hand side of Eq. ( pT| ) 
will counteract the flrst one. The existence of a nonvanishing bulk viscous pressure implies that the temperature 
decreases less than for a perfect fluid, tt tends to 'reheat' the cosmic matter. 

Use of the equations of state ( [l(]| ) and the temperature law (jlj) allows us to write the evolution equation (||) for 
the viscous pressure as 
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is the square of the adiabatic sound velocity Cg. The quantity Cf, is the propagation velocity of viscous pulses, 
determined by [p|j4l| 
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(14) 



The sound in a viscous medium propagates with a subluminal velocity v where pl| 

v'=cl + cl<l. (15) 

With the help of the field equations {k is Einstein's gravitational constant and H is the Hubble parameter H = 0/3) 



Kp = 3H\ H = --{p+p + tt) 



(16) 



for an homogeneous and isotropic, spatially flat universe, we may write 



KTT = -3-fH^ - 2H 



where 7 = 1+ p/p, and 
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Combining the last two relations with Eq. (12), the viscous pressure and its derivative may be eliminated to yield 
the basic causal evolution equation for the Hubble parameter: 
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i7+-7ij2 = o. (19) 



This equation is valid for arbitrary, even arbitrary time varying equations of state. To get an idea of the physical 
implications of Eq. (|l9|) we introduce the notations 



P A - , 1 , "^^ n - , dp dp 2 

Td^Jdf^^^ = "''' + '+d-p^''-^'-'-d-p^^-d-p-'^^ 
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with 1 < 7 < 4/3 and assume the values of r, A, B, C, and Cf, to be constant. In the ultrarelativistic limit with 



7 = 4/3, dp/dp 



1/3 we have r = I, A = 2, and B — C — 0. For non-relativistic matter, characterized by 



7 = 1, dp/dp = 2/3, c^ <C 1 the corresponding quantities are r ^ 1, A ^ 1, i? ::^ 1, and C ~ 2/3. The general 
equation (|l9[) then simplifies to 
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This equation has stationary solutions H = for which 
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is valid. Causality restricts Cb to values c^ < 1 



' [cf. Eq. (|l5|)]. For equations of state close to that for non- 
relativistic matter one has c^ < 1 as well as _B 3> 1 and C ~ 2/3 and the relaxation time r is much shorter than 
the cosmological time scale H^^. However, for equations of state close to that for radiation with r « 1, c^ < 2/3, 
B « 0, and C « 0, the relaxation time may well be of the order of the Hubble time, i.e., the nonequilibrium may be 
"frozen in" . Such kind of "freezing in" may be regarded as a necessary condition for successful inflation |40[| since 
otherwise the system would relax to equilibrium in less than one Hubble time. Since, on the other hand, one expects 
T to be typically of the order of the mean free collision time and this collision time itself has to be smaller than the 
Hubble time for a fluid description to make sense, we face here serious applicability limits of a viscous fluid approach 
to inflationary cosmology. In fact, the situation is still more unfavorate. The stationary case implies ir ~ ~ {p + p) 
and T oc a3|7r-dp/dpj ^[^^i a ex exp [Ht] , i.e., the temperature is exponentially increasing for any equation of state. 
The number density in contrast is exponentially decreasing as n ex a^^. This decrease compensates the increase in 



T, resulting in a constant energy density p. A behavior such as this, although formally not excluded, seems highly 
unrealistic. 

Let us now look for solutions a ex i^ with constant values of g, equivalent to H = qt^^, H = —qt^^, and H = 2qt^^. 
The parameter tH becomes 

^H = ^jJ^Q , (23) 



with the abbreviation 



3ci 



^ cl + l{B-C) + ^[A-2^3q{B + C)] ■ 

In the ultrarelativistic limit Q approaches Q = 1, while Q ^ 1 in the opposite limit of non-relativistic matter. On 
the other hand, we find 

-^ = ^-^ ^ rH = -^^Q. (24) 

p + p 3jq 3cg p + p 

For TT 7^ we have r — > in the limit c^ ^ oo which is the Eckart case. The last equation of Eq. ( |24| ) represents an 
algebraic relation between the relaxation time r and the viscous pressure n. Solving the first relation of Eq. ( P4| ) for 
q yields 
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Eqs. (p3) and ([25|) imply a constant ratio 7r/p under the condition of a power-law behaviour of the scale factor. The 
temperature law (UW becomes 
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The case q = 2/(37) corresponds to tt = (perfect fluid limit). Since Q > only negative values of tt guarantee 
a positive relaxation time r which is consistent with a positive bulk viscosity coefhcient according to Eq. (|lj) and 
positive entropy production according to Eq. (O). The case q = \ corresponding to a = 0, is characterized by 
■K = — (1 — 2/ (37)) {p + p). The energy density decays as p oc a^^, i.e., it mimicks the behaviour of so-called "K- 
matter" Q]. The temperature increases, however, as T ex a^''''''^^"^/'^'^'''-'-'"^''/'^''! in this case. The number density 
n decays as n ex a~^ for arbitrary q. For any g > 1 we have accelerated expansion, i.e., power-law inflation. For 
very large q the viscous pressure n approaches — {p + p), the case of exponential inflation. For any value q > 2/3 the 
non-equilibrium pressure exceeds the equilibrium one, i.e. |7r| > p, and we are beyond the strictly justified range of 
applicability of the theory. 

The considerations so far summarize and generalize previous presentations of bulk viscous cosmology pq]. They 
confirm that within the framework presented above, relying on a constant number of particles, the concept of bulk- 
viscosity-driven inflation is hardly convincing. However, particle number conservation is an exceptional situation in 
early cosmology. In the following section we shall extend the theoretical basis of our investigations insofar as we will 
admit the possibility of a nonvanishing source term on the right-hand side of Eq. (0) , i.e. , particle number conservation 
will no longer be required. 

III. THE UNIVERSE WITH PARTICLE PRODUCTION 
A. General theory 



In case the fluid particle number is not preserved, the conservation law (pi) has to be replaced by the balance [17 



iV.'^ = ri + en = nF , (27) 

where F = N/N is the change rate of the number N = na^ of particles in a comoving volume a"^. For F > we have 
particle creation, for F < particles are annihilated. A nonvanishing particle production rate is known to generate 



an effective bulk pressure of the cosmic fluid. This is most simply demonstrated for "isentropic" particle production. 
From the Gibbs equation 

Tds^d--pd- (28) 

n n 

together with the balances (0) and (p^) we obtain 

nTs = -Gtt- (p + p)r . (29) 

"Isentropic" particle production is characterized by s = 0. Under this condition the equilibrium entropy per particle 
does not change as it does in dissipative processes of the type discussed in Sec. II. Instead, the viscous pressure is 
entirely determined by the particle production rate |17| , |l8| : 

s = ^ n^-{p + p)^. (30) 

The cosmic substratum is not a conventional dissipative fluid but a perfect fluid with varying particle number. There 
is, however, entropy production due to the enlargement of the phase space of the system. First we consider the 
analogue of the Eckart theory (subscript E). From the entropy flow vector S'^ = nsu"" one obtains a nonvanishing 
entropy production although s = 0: 

S%a = nsV . (31) 

While in principle the rate F has to be calculated from microphysics (it may also be prescribed by symmetry require- 
ments connected with specific inelastic self- interactions of the cosmic matter [ plp^ j ) , a phenomenological determina- 
tion along the lines of linear, irreversible thermodynamics is possible here as well. To this purpose we rewrite Eq. 
( pi| ) with condition (|3^) as 

^^■'^- p + pT- ^-^^^ 

(Notice that nsT ^ p + p for a vanishing chemical potential.) The simplest way to guarantee 5'^.^ > is the linear 
ansatz 

Ts = -C^e^-Ce//e =^ rE = ^e\ (33) 

p+p p+p 

with a coefficient ^ > 0. Of course, this implies F^ > 0. Only a positive Fg is compatible with the second law of 
thermodynamics. The existence of relations of the type (^3|) shows that there is a close analogy to the conventional 
first-order theory of dissipative processes (in which F = and s j^ 0). It appears natural to extend this analogy also 
to the second-order theory where the viscous pressure becomes a dynamical degree of freedom. In the present setting 
this allows us to dynamize the particle production rate F which was introduced as an additional parametern in the 
balance (27). Our approach starts again with the expression (^ for the entropy flow vector. However, due to the 
possibility of a varying particle number, Eq. (ph is now replaced by 

5™ =ni + s7V™-l(^^^.™) tt'-^ttt:. (34) 

While the expression (ph was derived under the condition of small deviations from the standard perfect fluid behavior, 
we will not explicitly apply this restriction in the following. Instead, we will show that in case n is entirely a 
consequence of cosmological particle creation, the dynamics based on Eq. (p4) may provide its own limits. Such kind 
of self-limitation does not exist for a conventional viscous pressure tt. Nevertheless, a really satisfactory theory should 
be based on a higher-order generalization of the expression (||) . (For a recent attempt in this direction see |Q ) . The 
present approach may be regarded as a flrst step which exploits the fact that the entropy flow (^ is the simplest 
option which allows us to dynamize the particle creation rate on the phenomenological level at all. 

Focusing on the case (|30| ) of "isentropic" particle production, i.e., on processes with constant entropy per particle, 
means that the variable tt does not characterize a conventional nonequilibrium but a state with equilibrium properties 
as well, which, however, deviates from the specific equilibrium with F = 0. The equilibrium condition s = is 
maintained through the entire stage for which F > 0. In other words, the system may always be regarded as a perfect 
fluid [see the comment following Eq. (po|)]. This implies that the variable r, which plays the role of a relaxation time 



here as well, has to be interpreted now in a different manner. While the quantity r of the previous section (Sec. II) 
describes the relaxation from states with s > to those with s = due to internal (elastic) collision processes and is 
of the order of the mean free collision time, we have now a "relaxation" from states with s — and F > to states 
with s = and F = 0. The corresponding process is not related to interactions in which the mean free collision time 
plays a role. Consequently, it will no longer be contradicting to have a relaxation time of the order of the Hubble 
time and, on the other hand, internal collision time scales which are small enough for a (perfect) fluid description to 
be sensible. This picture relies on a separation of interactions which guarantee a perfect fluid behavior from those 
which are responsible for the particle creation. The time scales of interest here refer to the latter, i.e., to the particle 
number non-preserving part of the interactions. The time scales of the former are assumed to be much smaller and 
will not explicitly appear in the following. With this basic difference in mind we will continue to use the symbol r to 
characterize the "relaxation" from F > to F = 0. 

Under the condition of "isentropic" particle production, the change rates for the relevant thermodynamic variables 
are 
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and 



p = -{e-T)ip+p) , p = -cUe-T){p + p) . 

With the help of these relations the entropy production density (p^) may be written as 
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Inserting here for the last term 
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which follows with the help of Eqs. (H), (pq), and (pq), the expression (pq) becomes 
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Analogously to the standard procedure of Sec. II, the simplest way to guarantee 5.™ > is again a generalized linear 
relation 
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with C > 0, which implies 
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for the entropy production density. Use of Eq. dlj) in the generalized linear relation (BO) allows us to obtain a 
dynamical equation for T/Q: 
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which is the counterpart of Eq. ( ]12| ) of Sec. II. For r ^- one recovers the Eckart-type relations 
The field equations (16) and Eq. ( |30| ) may be combined to yield 



r _ 2 H 

Use of the last relation in Eq. (p2[) provides us with the following second-order equation for H, 
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which is the counterpart of Eq. (|T^ of Sec. II. For p = p/3 and ns = {p + p) /T the general equation 
reduces to 



(45) 
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which is the analogue of Eq. (pT|). Again there exist stationary solutions H — 0. From the relationship (44) we have 
F = 3H in this case, equivalent to constant values of n, T, p and p according to the evolution laws (|3^) and (|36|). The 
apparently unphysical behaviour of the temperature of the previous section (see the discussion following Eq. (p^)) 
docs no longer occur. Since n is stationary as well, there is no exponential dilution of the fluid either. This behavior 
demonstrates that the extension of the fluid dynamics to not necessarily small values of tt is more satisfactory in case 
IT represents cosmological particle production. For the parameter tH we find 
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Comparing the last expression with Eq. ( p^ ) for i? = C = shows that particle production leads to a higher value for 
tH than "conventional" causal cosmology, i.e., the relaxation time is enlarged. (We recall that the "relaxation" time 
T here refers to a different process than in Sec. II). On the other hand, the applicability of our formalism is restricted 



to ct > 



The general causality restriction 



< 1 



< I [cf. Eq. (O)] is compatible with c^ > i. For values 



c^ < 2 there does not exist a physically meaningful stationary solution of Eq 



the relativistic matter dynamics for different cases on the basis of Eq. (46) 



In the following we shall study 



B. The linear case T <x H 



Looking for power-law solutions a ex i^ of Eq. (46) one obtains 
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The particle production rate becomes 



r=(i_-)3i. 



(48) 



(49) 



i.e., T (X H follows as a consequence of the power-law ansatz, while the expression (|25|) remains valid with tt given by 
Eq. (pO). This allows us to write the parameter tH as 
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or, in terms of F, 



tH ^ -— 7. (51) 
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For r — > 37?, corresponding to g ^ cxd we recover the case ( [471 ) . Comparison of the parameters ( pOf ) and (23) for 
7 = 4/3 and (5 = 1 shows that also in the present case the relaxation is generally larger than for conserved particle 
numbers. Since — F is 

e-F^Aif, (52) 
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the balances ( |35| ) and ( |3q ) reduce to 
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resulting in 
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n lyi a ^1 , T cc a ^1 , p cc a 1 . (54) 



For q = ■^ we recover the standard perfect fluid behaviour. For q > ^, equivalent to F > all quantities decrease at 
lower rates than ior q — ^. Unlike the F = case of Sec. II, the equilibrium relation n (x T^ continues to hold here 
for any value of q. An essential advantage of the present setting is an inherent self-limitation of the theory, connected 

with the denominator cl — ^ [^) ~ ^l ~ "^ \~) ™ E^- (0)- ^Jur considerations make only sense for ( - j < ^c\. 
Otherwise we would obtain a negative r which is clearly unphysical. This represents an example for the property 
mentioned below Eq. (p4|), according to which the magnitude of the effective viscous pressure may be limited by the 
internal fluid dynamics. 

C. The quadratic case F oc H^ 

To study more general situations it is convenient to replace H in Eq. ( |44| ) by 

H = -—a = H'Ha , 
da 

where H' = AH /da. The resulting equation is 
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For F — aiH with a — const, it integrates to 
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and 
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a ex i^5MT^^ (F = aSH) , (57) 

which are the previous results for the linear case: A linear relation between F and H implies a power-law behaviour 
of the scale factor and vice versa. A more interesting case is a dependence F ex H^ , which was already discussed in 
p9[ in the context of matter creation due to a decaying cosmological constant: 

T H 

with P = const. The quantity He = H {a^) is the Hubble parameter at some flxed epoch with a = Og. With the 



quadratic ansatz (58) integration of the differential equation (|55|) yields 



H 



-iJ.. . 



{l~P)ah +pal 



(59) 



For a — !■ we have H —> P ^H^ = const, i.e., an exponentially accelerated expansion (a > 0), while for a ^ Og we 
recover the familiar dependence H oc a^^^/^ for FLRW universes (with d < 0). If we identify the intermediate value 
of a at which d = with Oe, we have He = —H^. From Eqs. (E3) and (BSt) it follows that (3 is then fixed according to 



i.e., for relativistic matter we find 



/3 



37 



H 



at 



tH, 



a^ + a% 
While a linear dependence of F on H corresponds to a power-law solution for the scale factor, we have now 

1 



t = te + 



AH, 



In 



1 



(60) 



(61) 



(62) 



This is the same dependence as the one obtained on the basis of a phenomenological model for a decaying cosmological 
constant by Gunzig et al. [ p9| . (For alternative thermodynamic descriptions of cosmological vacuum decay see p6| , p3| ) . 
The resulting cosmological scenario, however, is different here from the model in |29|] . In our case F is limited to 
r < 3iJ, while in [g^ the creation rate may be larger than this value. This gives rise to a difference in the thermal 
evolution at very early stages (see below). But our main point here is the following. In obtaining the solution (|6l| ) 
we have not really solved our basic equation (Esh but only Eq. (E3) under the assumption (pq). The remaining part 
of the dynamics is given by Eq. (|4^). In other words, we are dealing with an enlarged set of basic equations. For 
p = p/3 and vanishing chemical potential Eq. (B3) may be written as 



Qt 





r 
e 




c^ - 


-H^r- 


(r/e)- 
" e 



(63) 



It is obvious that the dynamical degree of freedom F/0 determines the parameter Qt as well. Since we prescribed a 
specific functional dependence for T /H, Eq. ( |6^ ) fixes the parameter tH. On the other hand, in order to be physically 
meaningful (positive relaxation time), any choice of F/O is restricted to 



(r/e)- 
e 



< ct 



(64) 



This condition on F (or, equivalently, on the corresponding effective viscous pressure according to Eq. (pO[)) is a new 
result. In the present case the left-hand side of the last relation becomes 



1/ry (F/ey 
2 \e) e 

It follows that the parameter tH is given by 

tH 



The condition (p4) specifies to 



H < 



1 H 

~iH~e 



1_H_ 
8 He 



1 ] 






H 


6 „2 , 1 if 
^b ^ 3 ffe 


1- 


7 H 

8 He 


He 



(65) 



(66) 



21 



He 



(67) 



This represents a restriction of the admissible Hubble rate in dependence of the parameter c^. Since H starts at 



Hq = 2He [cf. Eq. mm for 



0], this inequality is always satisfied for 



> ct > 



i.e., our whole scenario is 
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compatible with general causality requirements for effective bulk viscous fluids from the beginning. For the maximally 
acceptable value of c^ which is compatible with the causality requirement, c^ = 2/3, (see the discussion following Eq. 
(p7|)), the initial value of r is 



To 



2^0-' 



(68) 



i.e., the relaxation time exceeds the Hubble expansion time. For c^ < i the theory is only applicable at later 
periods of the evolution but not at very early times. Even for c^ <C 1, however, where condition (jG^) becomes 
H < jHe, the applicability still starts within the phase of accelerated expansion {H > H^). These restrictions are 
again manifestations of the fact that the theory considered here provides its own limits of validity. At the same time 
this precises the conditions under which a fluid picture of the inflationary universe is permitted. 

The results obtained so far give rise to the following conclusions. As long as one assumes a power-law solution for 
the scale factor, F is necessarily proportional to iJ, i.e., the ratio T/H is constant. The impact of a nonvanishing 
particle production rate on the cosmological dynamics does not change in time under this condition. Given F is large 
enough at a certain time to generate a power 9 > 1, i.e., accelerated expansion, there is no way to terminate such a 
phase. Consequently, models with a ex t' which imply T/H = const may be realistic at most piecewise. For a model 
with F ex H^, on the other hand, the time dependence of the scale factor is more involved (see Eq. (^)) but there 
is a variation in the ratio T / H . The impact of a nonvanishing particle production rate on the cosmological dynamics 
is different in different epochs. In particular, it implies a transition from an inflationary phase to a conventional, i.e. 
decelerating FLRW period. This transition is accompanied by a decrease of the parameter tH . For radiation the 
latter quantity is generally given by Eq. (^6|) . For H ^ Hf,, equivalent to a ^ Ce the ratio tH tends to zero and the 
cosmic substratum approaches a perfect rclativistic fluid behavior with T/H <C 1, i.e., negligible particle production. 

Combination of the Friedmann equation (16) with the Hubble parameter ( |6l[ ) allows us to obtain p — p (a). 



Sifg 2 



(69) 



where we have replaced k by the Planck mass rup according to k = ^tt /mi. In the following we shall assume | > c^ > i 
to ensure that condition ( |67|) is always satisfled. Initially, i.e., for a ^ 0, corresponding to i — > —00, the energy density 
is constant, while for a^ a^ the familiar behaviour p oc a"** is recovered. In our approach the evolution of the universe 
starts in a quasistationary state with a finite initial energy density, while the two-component model of Gunzig et al. 
p9[ implies a vanishing radiative energy density at a = for i — > — cx3. This apparent contradiction is easily resolved, 
however, if we split the energy density (|69|) according to p = pi -\- p2 with 



Pi 



ml 

27r 



P2 






-\3 



The part p2 corresponds to a cosmological "constant" in ||29|| which decays as a"® for a ^ Oe , while the part pi describes 
rclativistic matter with pi — > for a -^ 0. While this makes obvious that an interpretation of our approach in terms 
of a decaying cosmological "constant" is also admitted, the dynamics of particle number density and temperature is 
different from [E9[ . In our case these quantities behave as 



and 



n = uq 



T = To 




2V2r. 



\/2Te 



-,3/2 



-,1/2 



(70) 



(71) 



respectively, where the subscripts and e refer to the corresponding values at a = and at a — Ue- All the 
thermodynamic variables remain finite for a —f and they approach the correct dependences for radiation for a 3> flg. 
In particular, the universe starts with a finite maximal temperature, while the initial radiation temperature in |29|] 
is zero. In a sense, the radiative matter itself plays the role of a decaying cosmological "constant" in our approach. 
During the interval —00 < t < te, equivalent to < a < Oe, the energy density decreases by a factor of 4, while n and 



T decrease by factors of 2-^2 and \/2, respectively. For the number N 
find 



of particles in a comoving volume a we 
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N = 2V2Ne 



3/2 



3/2 



N 



f 



(72) 



where N^ — nf,a\. The particle number A'' starts from TVq = at a = and subsequently grows to the final constant 
value Nf — 2^/2Ne for a ^ a,,. This value may be identified with the number of particles in the presently observed 
universe which is of the order 10*^. Since the epoch a = ae with a (te) = the total particle number has grown by a 
factor of 2\/2- The circumstance that particles are absent initially (A^o = 0) is consistent with the interpretation of 
the corresponding initial energy density 



Po 



3H, 

2tt 



-Hi 



(73) 



following from Eq. (B9) in the limit a — > 0, as a vacuum energy density. Consequently, our approach may be regarded 
as equivalent to the model of a decaying vacuum. 

We finish this section by considering the behavior of the entropy in our model universe. With Eqs. (M), (pffl 
and (plj), the entropy flow vector (|5|) may be written as 



5-^ = 



(74) 



where s* is 



2c2 



I H 
2I1, 



2c2 



(75) 



This quantity is an analogue to the non-equilibrium entropy per particle for viscous fluids with conserved particle 
number. In the present context it is more adequately called an effective entropy per particle. For a — > we have 
(recall that s = 4) 



1 

24 



(76) 



This 
Eq. ( 



mantity is positive for c^ > 1/2 which coincides with the condition for the relaxation time to be positive [cf. 



a) for H = Ho = 2He ]■ 



With cl 



2/3, the maximally allowable value, we get 



cl 



(77) 



The effective entropy s* per particle is reduced considerably compared with the value s = 4 for F = . A reduction to 
s* = 1 is compatible with the causality requirement. During the period where the particle production is highest, i.e., 
for H > He, the effective "nonequilibrium" contribution to the entropy is highest as well. For H ^ He corresponding 
to a 3> fle, the quantity s* approaches s. While s is constant according to condition (30), s* is not. From the definition 
([Tq) we obtain 



16 a^ 

cl a^ + al 



.2 i3 



a^ + ai 



He 



(78) 



This expression grows as s* ex a? for a <^ ae and decays as s* ex a ^ for a 3> fle. The effective entropy in a comoving 
volume is given by 



S* = ns*a-^ = Ns* 



(79) 



It increases both due to the increase in A'' from Aq = to Nf and due to the increase of s* to s*^ — s — A for a 3> flg. 
We emphasize that different from "standard" inflationary models which rely on the dynamics of a scalar field, this 
scenario does not imply a separate reheating phase. The mechanism which drives inflation is at the same time entropy 
producing. In the following Section we establish a relation between the fluid model of the universe presented so far 
and models which are based on the dynamics of scalar fields. 
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IV. SCALAR FIELD DYNAMICS 



Inflationary cosmology is usually discussed in terms of scalar fields with suitable potentials. The usual picture 
consists of an initial "slow roll" phase with a dynamically dominating (approximately constant) potential term, which 
generates a de Sitter like exponential expansion, connected with an "adiabatic supercooling" . During the subsequent 
"reheating" phase which in itself is of a highly complicated non-equilibrium structure, the scalar field is assumed to 
decay into "conventional" matter. The entire entropy in the presently observable universe is produced during the 
reheating phase according to these scenarios. As to the mechanism of entropy production the "standard" (i.e. based 
on scalar field dynamics) infiationary picture essentially differs from the fluid particle production scenario of the 
present paper, since here the entropy is produced already during the period of accelerated expansion. On the other 
hand, as far as the dynamics of the scale factor is concerned, there exists a close correspondence between the role of an 
effective negative fluid pressure (due to the production of particles) and a suitable scalar field potential. What counts 
here is the magnitude of the effective negative pressure, independently of its origin (scalar field potential or particle 
production). What one would like to have is an interconnection between both these different lines of describing the 
early universe. This would allow us to switch from the scalar field to the fiuid picture and vice versa. It is the purpose 
of this Section to show how the effective imperfect fluid dynamics discussed so far in this paper may be translated into 
the dynamics of a minimally coupled scalar field. In particular, we shall find a scalar field equivalent for the previously 
considered [cf. Sec. Ill] comoving fiuid entropy which was shown to be increasing in time due to cosmological particle 
creation. To this purpose we start with the familiar identifications 



P=\^'' + V (0) 



P 



P = p + Tr=^^^-V{^) , 



p-P = 2V{<l>) 



(80) 
(81) 



With the expressions (p!l, (M), (pq), and (61), we obtain for p + P 



p + P = 4^—-m~ 



a^ 



at 



(82) 



For (f> we write = (f>'a — (p'aH, where 0' = dcj)/da. It follows that 

rriri 1 



Integration of this equation yields 



y = ±- 



bp 



2n ^a2 + al 



(83) 



f>o = 



rup a a . , 

lArsh — =^ — = smh 



\/27r fle a. 

where we have restricted ourselves io 4>> 4>q. The potential is given by 



27r- 



TOd 



(84) 



P-P Hi ^ 

= TTit, 

2 27r P 



+ 3 



n3 



a^ + at 



V , 



(85) 



0-2 cosh^ 
y _ ^e_2 



-to: 



2^^ ^ cosh^ 



/2^i 



2 tt2 3-2 tanh^ 

-"e 2 

— = — -nio 

2^^ cosh^ 



/27r^ 



27r^ 



2H, 



(86) 



This expression was found by Maartens et al. p7| as exactly that potential which produces the dynamical behaviour 
(p2) of the scale factor. The rate F as a function of (j) is 



QHe 



cosh 



(87) 



In terms of the scalar field the cosmic evolution starts with 



13 



Vo 



3^e 2 



and 00 = 0. At a = Cg the value of the potential is reduced to 



47r ^ 



(89) 



The relations found so far in this section allow us to change from the fluid to the scalar field picture and vice versa at 
any time of the cosmological evolution. In particular, they may be used to express the entropy d79) in terms of the 
scalar field. Namely, combining Eqs. (|7^), ( [75|) and (|7|) with the second relation of Eq. (84), we obtain 



S* = ANf tanh' 



1 7T cosh 

4 



'2-K- 



(90) 



as the "effective" scalar field entropy of a comoving volume. This amount of entropy may be attributed to the universe, 
independently of the specific matter model. The corresponding first-order limit S for c^ ^ oo is 



S* 



S = ANf tanh 



^2n 



0-00 



(91) 



The property of the scalar field configuration to represent a certain amount of entropy may be used to avoid solutions 
of the type of "jump" solutions p7| according to which an instantaneous reheating due to a discontinuity in the 
effective pressure at the transition from an early scalar field epoch to a subsequent radiation dominated FLRW phase 
is supposed to generate the entropy in the universe. 

We mention that also for a power-law behaviour a oc f^, equivalent to the linear dependence ([49|), an equivalent 
scalar field dynamics may be formulated. By similar steps like those used to derive the expression (bq), we find for 
the potential in the linear case 



V 



(3g- 



-T^exp 



tf 



1 



1 
Yq 



3H 



(92) 



where ti is some initial reference time. The potentials (p6[) and ( |92| ) are of the type of exponential potentials. 
Exponential potentials arise naturally in particle physics [ |48| , [49| and are known to be typical for power-law solutions 
for the scale factor [p0|-p3[ . They were also studied in the context of interacting scalar field cosmologies 1 5J| and in 
discussions of the recent SNIa results according to which the present universe might be accelerating |^,^ | . 

The option to use either the fluid or the scalar field picture and to switch between them gives rise to a complementary, 
more comprehensive overall picture of the dynamics of the universe. Moreover, it may contribute to establish a link 
between fluid cosmology and particle-physics-motivated investigations of the early universe. 

As a final remark we point out that rewriting the fluid dynamics into the dynamics of a scalar field implies only 
those subset of equations which does not depend on the parameter tH. There is no counterpart of relation (pq) in 
"conventional" scalar field theory. Compared with the latter, the causal fluid dynamics is more restrictive. It provides 
us with limits for the admissible inflationary dynamics which do not exist for scalar fields. Again, this may be traced 
back to the inherent self-limitation of the causal theory for isentropic particle production. 



V. CONCLUSION 

We have developed a setting according to which isentropic cosmological particle production is a dynamical process 
described with the help of the second-order Israel-Stewart theory for (effectively) imperfect fluids. This picture uses 
the well-known equivalence between an effective viscous fluid pressure and cosmological particle creation to make the 
rate for the latter a dynamical degree of freedom. The resulting dynamics was shown to be inherently self-limiting, 
i.e., it may restrict the magnitude of the effective viscous pressure. We have precised the conditions under which an 
effective bulk pressure, compatible with the general causality requirements of second-order thermodynamics of the 
Israel- Stewart type, may drive a phase of inflationary expansion. In particular, for a quadratic dependence of the 
particle production rate on the Hubble parameter we have established a scenario according to which the evolution of 
the universe starts in a de Sitter phase with finite, stationary values of density and temperature, and subsequently 
evolves smoothly to a standard FLRW stage. The corresponding relaxation times turned out to be of the order of the 
Hubble time during the phase of accelerated expansion. This scenario was also shown to be consistent with the model 
of a vacuum energy which decays into relativistic matter. Furthermore, we have clarified that part of the effective 
imperfect fiuid dynamics is equivalent to the dynamics of a minimally coupled scalar field. This correspondence was 
used to find a scalar field equivalent for the (generally increasing) comoving fluid entropy. 
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